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A particular example is produced to prove that quantum walks can be used to simulate full- 
fledged discrete gauge theories. A new family of 2D walks is introduced and its continuous limit is 
shown to coincide with the dynamics of a Dirac fermion coupled to arbitrary electromagnetic fields. 
The electromagnetic interpretation is extended beyond the continuous limit by proving that these 
DTQWs exhibit an exact discrete local U{1) gauge invariance and possess a discrete gauge-invariant 
conserved current. A discrete gauge-invariant electromagnetic field is also constructed and that field 
is coupled to the conserved current by a discrete generalization of Maxwell equations. The dynamics 
of the DTQWs under crossed electric and magnetic fields is finally explored outside the continuous 
limit by numerical simulations. Bloch oscillations and the so-called E x B drift are recovered in the 
weak-field limit. Localization is observed for some values of the gauge fields. 
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I. INTRODUCTION 

Discrete Time Quantum Walks (DTQWs) are formal 
generalizations of classical random walks. They were first 
studied in a systematic fashion by D. A. Meyer [1], while 
similar but different quantum discrete dynamics was first 
considered in [mi- DTQWs have been realized experi¬ 
mentally with a wide range of physical objects and setups 
EHn, and are studied in a large variety of contexts, 
ranging from fundamental quantum physics El ng to 
quantum algorithmics [Hin], solid state physics [T5HIH] 
and biophysics [191EO] • 

Particular quantum cellular automata EH HD , among 
which DTQWs [23] defined on various regular lattices 
EH, are known to reproduce, in the continuous limit, 
the dynamics of free Dirac fermions in one, two or three 
spatial dimensions. Recently, such connexions have been 
extensively extended to Dirac fermions coupled to gauge 
fields [25ti3^ . More precisely, ID DTQWs have been pro¬ 
posed which reproduce the dynamics of Dirac fermions 
coupled to arbitrary electric [25H27] and/or gravitational 
[28H3T] fields, and a 2D DTQW simulating the coupling 
of a Dirac fermion to a constant uniform magnetic field 
has been proposed in [32] . 

For all existing DTQWs, the gauge fields are encoded 
in the time and space dependence of the operator ad¬ 
vancing the fermion in discrete spacetime. They act on 
the fermion but the dynamics of the fermion has no ef¬ 
fect on the gauge fields. In other words, the gauge fields 
play the roles of imposed external fields. In particular, 
they are not advanced by their own discrete dynamical 
equations, as is for example the case in Lattice Gauge 
Theories (LGTs). The main purpose of this article is 
to remedy this problem and introduce the first complete 
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self-consistent model based on DTQWs where both the 
fermions and the fields are advanced by compatible dis¬ 
crete dynamical equations. 

The first brick in such a self-consistent model is a fam¬ 
ily of DTQWs which exhibit an exact discrete gauge in¬ 
variance associated to a certain group G and which de¬ 
scribes the coupling of Dirac fermions to arbitrary G- 
gauge fields. The existing literature contains only one 
such family. The associated gauge group is U{1) and 
the DTQWs describe the coupling of ID Dirac fermions 
to arbitrary electric fields. Electromagnetism is how¬ 
ever degenerate in ID. There is no magnetic field and 
Maxwell equations reduce to the Maxwell-Gauss equa¬ 
tion, which contains no time derivative. We therefore 
switch to 2D and introduce a new family of DTQWs 
whose continuous limit coincides with the dynamics of 
a Dirac fermion coupled to 2D arbitrary electromag¬ 
netic fields. We then show that these DTQWs admit 
(i) an exact discrete U{1) gauge invariance (ii) a gauge- 
invariant discrete electromagnetic tensor (f.e. gauge in¬ 
variant electric and magnetic fields defined on the dis¬ 
crete lattice of the DTQWs) (iii) a discrete conserved 
current. We finally combine the discrete electromagnetic 
tensor and the discrete conserved current into discrete 
gauge-invariant Maxwell equations which imply current 
conservation. This literal material is complemented by 
numerical computations which explore how the DTQWs 
which serve as a basis for the whole construct behave 
outside the continuous limit. Even outside this limit, the 
DTQWs display in the weak field regime several well- 
known features usually associated to standard continuous 
motions in electromagnetic fields, including Bloch oscil¬ 
lations and the so-called E x B drift. In the regime of 
strong fields, the discrete dynamics depends crucially on 
whether the fields are rational or not, as expected from 
previous work on other DTQWs [26l l27l [33] . We finally 
discuss possible applications of our results to quantum 
simulation and quantum algorithmics and highlight how 
the new discrete gauge theory based on DTQWs differs 
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from standard Lattice Gauge Theories (LGTs). 


II. ELECTROMAGNETIC DTQWS 

A. The walks and their formal continuous limit 

We consider DTQWs with two-component wavefunc- 
tions {2D coin or spin space) defined on a discrete (1 + 2)- 
dimensional spacetime where instants are labeled by the 
index j G N and space points on the 2D square lattice 
are labeled by the indices (p, g) G The evolution 

equation for the wave-function reads 

= U(6> (e^, m), T 2 (1) 


X U(0+(e„,TO),e^7,p,q 

, 0 ) T; 

1 ’ 


where the action of the shift operators Ti 

and T 2 

on 

the 2D wave-function 


)^ (the 

su- 

perscript T denotes the transposition) 

is: 






(2) 
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The coin operator U(^,^,(a) G U{2) is the product of 
three simpler operators: 


U((9, a) = D^lx C{0) X S(0 (3) 



0 


COS 0 i sin 0 


" o' 

0 

gia 


i sin 0 cos 0 


. 0 . 


The first operator S(^) is a spin-dependent phase shift 
parametrized by the angle the second operator C{0) 
is a standard coin operator with angle 0 and the third 
operator performs a global multiplication by the phase a. 

In the continuous limit, the parameter m which en¬ 
ters the definition of the constant angles 0^{e^^m) = 
if ~ interpreted as the mass of the walk 

and the three angles A^, A^, A^, which may depend on 
{jiPi q)^ will be interpreted as the components of an elec¬ 
tromagnetic potential. The positive parameters and 
eA are introduced to trace the importance of m and A 
and will tend to zero in the continuous limit. All the 
parameters and angles are dimensionless. 

The formal continuous limit of the DTQWs can be 
determined by the method used in [28l [29l [32j |34H36] : we 
first introduce a (dimensionless) spacetime-lattice step 
ei and interpret any (jf,p, g)-dependent quantity Qj,p,q 
as the value taken by a function Q{X^,X^,X‘^) at time 
Xj = jei and spatial position (Xf = pei^ X^ = qei). 
We then consider the scaling = e/ = e and let 

e tend to zero. We expand Eq. 0 at first order in e 
around the generic spacetime point (Xj^,X^,Xg). For 
the continuous limit to exist, the zeroth-order terms of 
the expansion must balance each other; this contraint is 
automatically verified by the DTQWs defined by 0 . 


The first-order terms of the expansion in e deliver the 
differential equation which determines the dynamics of 
the walker in the continuous limit. This equation reads: 


{i^^V^ — m)T = 0 . (4) 

Here, the 7 matrices are defined by 7 ^ = cri, 7 ^ = icr 2 , 
7 ^ = ids, where the Pauli matrices are 
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and = dp — iAp is the covariant derivative associated 
to Maxwell electromagnetism, with Aq = A^, Ai = — A^, 
A 2 = — A^. The 7 matrices satisfy the (l+2)-dimensional 
fiat-spacetime Clifford algebra 7 ^ 7 ^ + 7 ^ 7 ^ = 2r]^^l, 
where [ 7 ^^] = diag(l, -1,-1) is the Minkowskian met¬ 
ric. Equation 0 is the Dirac equation describing the 
dynamics of a spin 1/2 fermion of mass m and charge — 1 
coupled to the electromagnetic potential A. To consider a 
generic charge g, just perform the substitution A ^ —gA. 
In Eq. 0 - the characteristic speed is 1 because we have 
chosen the same value ei for the dimensionless time and 
space steps. 


B. Rate of convergence towards the continuous 
limit 

Since the formal continuous limit is obtained from the 
DTQW dynamics by keeping terms that are first order 
in e, one expects the discrepancy between a solution of 
the DTQW and the corresponding solution of the Dirac 
equation to scale as e^. This can be tested numeri¬ 
cally by computing the distance between an exact time- 
independent solution of the Dirac equation and the time- 
evolution of this solution by the DTQW. For simplicity 
we choose Aq = —EX^.Ai = 0 , A 2 = —BX^, which does 
not depend on X^ and generates crossed, constant and 
uniform electric and magnetic fields. The Hamiltonian 
and the momentum in the direction can then be di¬ 
agonalized simultaneously. For E = 0, the eigenstates are 
called relativistic Landau levels [32] . For 0 < /3 = E/B < 
1 , the eigenstates can be obtained from the relativistic 
Landau levels by a boost of velocity (3. The resulting 
eigenstates 0 /,k(X1,X2) = ^i{X^,K) exp{iKX^) and 
eigen-energies £i^k are labelled by a couple (/,X) where 
/ = 0 or / = (±,n) with n G N* and K is the eigen- 
momentum in the X^ direction. In one time step e, 
the DTQW evolves 4>/(X^,X) into a certain function 
Wi{X^^K) which should be approximated, at first or¬ 
der in e, by W(X\X) = exp(-ifz,K x e)4>z(X\X). For 
each X, the distance between the two functions H/(', K) 
and Wi{',K), can be evaluated by 


Si{K) ^ 


||W(-,x)-w(.,x)|| 

\\Wl{;K)\\ 


(6) 
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where || • || stands for the norm of a position {X^- 
)dependent function ^ defined on the lattice: 


11^11 = 


Pmax (c) 

P— Pmax(c) 


( 7 ) 


where Pniax(e) scales as 1/e. Figures]^ and [^display how 
Si{K = 0) scales with e for various values of I = (+,n) 
and for various values of /3, having fixed B = 1 and m = 
1. These figures clearly confirm that 6i{K = 0) scales as 
e^ for a large range of e-values 



FIG. 1. Distance 6i{K = 0) as a function of e for various 
values of I = (H-,n). 



FIG. 2. Distance 6i{K = 0) as a function of e for various 
values of /3. 


III. DISCRETE GAUGE INVARIANCE AND 
ELECTROMAGNETIC FIELDS 

The discrete equations 0 are invariant, not only under 
a global phase-change of the spinor T, but also under the 
more general, local gauge transformation 

^ ^,P,, = ( 8 ) 

where the three ‘discrete-derivative’ (finite-difference) 
operators are defined by 


do — {L — di — Ai/ca, ^2 — A2Ei/eA, 

(9) 

with 

{d^Q)j,p,q — Qj-\-l,p,q 
{'^lQ)j,p,q = {Qj,p-\-l,q + Qj,p-l,q) /2 
{'^2Q)j,p,q = {Qj,p,q-\-l + Qj,p,q-l) /2 (10) 

{XiQ)j^p^q = (Qji',p+l,g ~ Qj,p—l,q )/2 
{A2Q)j^p^q = {Qj,p,q-\-l — Qj,p,q-l) /2 . 

This local gauge invariance is a discrete version of the 
standard continuous U{1) local gauge invariance asso¬ 
ciated to electromagnetism and displayed by the Dirac 
equation Q. A straightforward computation now shows 
that the three quantities Fqi, F 02 and F 12 defined by 

{Fpu)j,p,q = — {dyAp) ( 11 ) 

are gauge invariant. These are clearly discrete ver¬ 
sions of the usual electromagnetic tensor components 
Fpu = dpAjy — duAp. In particular, Fqi and F 02 represent 
respectively the two components Ei and E 2 of a 2D dis¬ 
crete electric field (parallel to the plan of the (p, g')-grid) 
and the component E 12 represents a discrete magnetic 
field Bs perpendicular to the plan of the (p, g)-grid. 


IV. GAUGE-INVARIANT CONSERVED 
CURRENT AND DISCRETE MAXWELL 
EQUATIONS 


Let 


^j,P,q = U(6»+(em,m),e^y p g,0) Ti (12) 

be the state of the walker after the shift along the p 
direction and the first coin operation. The spatial density 
associated to '^j,p,q is 


j,p,q 




j,p,q 




(13) 


Introducing notations p=|T p-l-|T+p and J = | T+ p — 
we can rewrite the previous equation as 


Pj,p,q ~ {^lp)j,p,q {XlJ') 


j,p,q 


(14) 
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The same computation carried out for = 

results in 

{Lp)j,p,q = {'^2^j,p,q — {^2J')j,p,q • ( 15 ) 

Inserting ( p^ into ( p!^ gives the discrete conservation 
equation 

i^p)j,P,q ~ {'^2'^lp)j,p,q ~ {'^2^1sJ)j,p,q ~ {^2sJ)j,p,q • 

(16) 

All operators defined in ( p!Q| ) commute with each other; 
in particular, S 2 A 1 = A 1 S 2 , so that the conservation 
equation can be written 

{^pJ^)j,p,q — 0 5 ( 1 "^) 

where the new finite-difference operators read, 

Do = do^ Di = di'E2i ^2 = A2/e^, (18) 

and the probability current on the square lattice is given 


jO 

j,p,q 

~ Pj,p,q ~ 

1 + 

'^ 3 ,p,q 

7I 

"^j,p,q 

~ *^j,p,q ~ 

1 ytp,, 1^ - 


/2 

j,p,q 

~ *^j,p,q ~ 

1 - 

^j^p^q 


In the continuous limit, the discrete conservation equa¬ 
tion 0 becomes the standard conservation equation 
dpj^ = 0 of the 2D Dirac current with 

Having identified the discrete current and the finite- 
difference operators involved in the discrete continuity 
equation 0 makes it possible to write the following 
simple discrete equivalent to Maxwell equations 

p,, , ( 20 ) 

which connects the discrete electromagnetic tensor 
{Fuv) j.p.q to the discrete current {J^)j,p,q- Indeed, 
Eq. ( [^ has the standard Maxwell equations as con¬ 
tinuous limit, and ensures the conservation of the dis¬ 
crete current {J^)j,p,q because it implies {DyJ^)j^p^q = 
{Dj^DpF^^)j^p^q^ which vanishes identically because op¬ 
erators Dp commute with each other and because 
{F^^)j^p^q is antisymmetric. 

V. SIMULATIONS OUTSIDE THE 
CONTINUOUS LIMIT 

We now focus on constant and uniform discrete electric 
and magnetic fields, for example E = Ewi and B = Bus 
where ui and U 3 are two unitary vectors respectively 
along the p- (or X^-)axis of the grid and perpendicular to 
the plane of the grid. A potential generating these fields 

is {^o)j,p,q ~ DpCi^ {^ 1 )j,p,q — 0? {^2)j,p,q — DpCi. 

Walks with B = 0 {resp. E = 0) will be referred to as 


E-walks (resp. H-walks). Walks with E 7 ^ 0 and B ^ 0 
will be referred to as EH-walks. 

Quantities of particular interest are the probability of 
presence of the walker 

I = p or q, its time-dependent /-mean {resp. /-spread), de¬ 
fined as the time-dependent average {resp. square-rooted 
average) value of / {resp. l^) computed with P as time- 
dependent probability law on (p, q). 

All computations are carried out with e^m = 1 , e/ = 1 
and the same simple initial condition: 'ip~{j = 0 ,p, g) = 1 
if (p, g) = ( 0 , 0 ) and 0 elsewhere, ' 0 ^(j = 0 ,p, g) =0 for 
all {p^q). The only remaining free parameters are e^E 
and caB. As will now be discussed, DTQWs for which 
both caE and e^E are much smaller than unity exhibit 
regimes which resemble continuous physics. DTQWs 
with larger values of e^E and e^E behave very differ¬ 
ently, and can even localize. 

Figure shows the time evolution of the p-mean for 
several E-walks. For e^E = 0, the p-mean varies linearly 
with time. This ballistic transport is typical of homo¬ 
geneous DTQWs, i.e. DTQWs whose coin operators do 
not depend on the spacetime point. Moreover, transport 
occurs towards negative values of p only because the ini¬ 
tial state has a vanishing -0+. For e^E 7 ^ 0, the p-mean 
oscillates in time around the value = —0.5 m with 
a period which coincides with the so-called Bloch period 
PbiocIi = 27r/(eAE) with an error smaller than one time 
step. Bloch oscillations were first predicted by F. Bloch 
[37] and C. Zener |38] for electrons moving in solids. They 
have been observed in 2D photonic lattices [39] and IE 
electric DTQWs [26l|40]. As e^E reaches a sizeable frac¬ 
tion of 27r, Teioch becomes of the order of a few time steps. 
Another oscillating mode with period of the order of one 
time step then appears and dominates the dynamics. 



j 


FIG. 3. Time evolution of the p-mean for E-walks with 
caE = 0 (black, solid), 0.02 (magenta, dashed), 0.04 (blue, 
dot-dot-dashed), 0.08 (green, dot-dashed) 0.16 (red, doted), 
0.64 (cyan, solid). The oscillating period is Tbioch = Stt/(caE) 
with an error less than one lattice site. 

Figure displays the probability densities at time 
j = 500 for several EE-walks with caB = 0.16. For 
caE = 0 (left), the walker is quasi-confined around the 
origin, with a typical radius which slowly increases with 
the time j and is, at each j, a decreasing function of caE 
(data not shown, see [32] for detail). When caE 7 ^ 0, the 
walker spreads in the g direction, up and down. The bot¬ 
tom front propagates with a speed which coincides with 
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EjB^ as supported by Fig. This corresponds to the 
classical so-called ‘E x B drift’ of a charged particle un¬ 
der crossed constant and uniform electric and magnetic 
fields (see, e.g., m)- The roughly circularly symmetric 
‘Landau profile’ obtained for caE = 0 seems to be trans¬ 
ported at the drift velocity. The behaviour of the top 
front is counter intuitive from the classical perspective. 
The top-front spreads with a speed which seems inde¬ 
pendent of caE. a very similar behaviour has already 
been pointed out in [42] for quantum particles moving 
under the influence of super-imposed electric and mag¬ 
netic fields in a 2D periodic potential with tight-binding. 



p 

FIG. 4. Probability density, at time j = 500, for FB-walks 
with caB = 0.16. From left to right, caE = 0, 0.01, 0.02, 0.03, 
0.04, and Pmax = 0.0943, 0.0578, 0.0209, 0.0181, 0.0178. The 
bottom front corresponds essentially to the classical E x B 
drift. 


spread decreases considerably. This weak F-field regime 
breaks down partially for = 1 and completely for 
caB = 7r/3, while the g-spreading is essentially enhanced 
for strong values of €aE. For €aE = 7r/2 and values of 
caB which are not rational multiples of 27r, the walk 
seems to be almost localized in q (this is also the case 
in p direction, data not shown). Fig. [^focuses on this 
apparent localization. In the long-time limit, the walker 
spreads ballistically for values of caB which are ratio¬ 
nal multiples of 27r. This ballistic spreading is consid¬ 
erably reduced (quasi-localisation) for caB = ir/A e 
and 7r/3 -h and the walk seems to really localise for 
CaB = 7r/2 + 5. The p-spread displays the same qualita¬ 
tive behaviours (data not shown). 


eaE 



FIG. 6. Evolution of the g-spread as a function of caE for 
FF-walks with magnetic field (i) caB = 0.16 at times j = 100 
(red, dot-dashed) and j = 500 (magenta, solid) (ii) caB = 1, 
at times j = 100 (black, dashed) and j = 500 (blue, dot-dot- 
dashed) (iii) F = 7r/3 ^ 1.047 (green, doted) at time j = 500. 



FIG. 5. Time evolution of the g-coordinate of the bottom- 
front local maximum of the probability density, for FF-walks 
with F = 0.16. From top to bottom, F = 0 (black), 0.01 (ma¬ 
genta), 0.02 (blue), 0.03 (green), 0.04 (red) and 0.05 (cyan). 
This maximum propagates in the direction of E x B (up to 
small oscillations in the p direction) and with speed F/F up 
to a 1% precision. 


Previous work on DTQWs coupled to electric or mag¬ 
netic fields 1211 EH Eg have shown that walks with field 
values which are rational multiples of 27r (‘rational fields’) 
follow very peculiar dynamics. Fig. displays the g- 
spread of FF-walks as a function of caE at two times 
and different values of caB. For caB = 0.16, which is 
not a rational multiple of 27r, there is a weak F-field 
regime (from = 0 to — 0.06) in which the g- 
spread increases essentially linearly with caE. This is 
the regime of Figures and For caE > 0.06, the g- 



FIG. 7. Time evolution of the g-spread for FF-walks with 
caE = 7r/2 and caB = 0.16 (black, solid), 7r/4 (red, dashed), 
7r/4 + £ (red, solid), 7r/3 (blue, dashed), 7r/3 + £ (blue, solid), 
7r/2 (green, dashed), tt/ 2 -\- e (green, solid), with e = 0.04. 


VI. CONCLUSION AND DISCUSSION 

We have introduced a new family of 2D DTQWs which 
coincides, in the continuous limit, with the dynamics 
of a Dirac fermion coupled to arbitrary electromagnetic 
fields. The wavefunction of these DTQWs has two com¬ 
ponents and the DTQWs explore the 2D square lattice 
by advancing alternately in each of the orthogonal di¬ 
rections. Similar, albeit simpler 2D DTQWs have been 
discussed for example in IMllla]- We have shown that 
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the new DTQWs introduced in this article possess an 
exact discrete local U{1) gauge invariance, a discrete 
gauge-invariant conserved current and a discrete gauge- 
invariant electromagnetic field, and that field and current 
can be coupled by discrete generalizations of Maxwell 
equations. We have also explored the behaviour of the 
DTQWs outside the continuous limit, under weak and 
strong fields. For weak fields, we have observed discrete 
versions of the Bloch oscillations and of the so-called 
E X B drift. We have also observed localization for some 
higher values of the fields. 

The results of this article prove that DTQWs can be 
used to build full-fledged discrete gauge theories and that 
laboratory experiments based on quantum walks can, at 
least in principle, simulate these theories (see for exam¬ 
ple [26] for a discussion of a quantum walk experiment al¬ 
ready carried out which simulates Dirac fermions coupled 
to ID electric fields). On the technical side, the construc¬ 
tion we have presented should naturally be extended, 
not only to Maxwell electromagnetism in AD spacetime, 
but also to other Yang-Mills gauge theories. Develop¬ 
ing second-quantized versions of these discrete theories 
should also prove interesting. 

A full comparison of possible discrete gauge theories 


based on DTQWs with the usual Lattice Gauge The¬ 
ories (LGTs) HH |45] is beyond the scope of this arti¬ 
cle. Let us simply mention two differences. First, unlike 
the ‘U’ parallel transporters in LGTs, gauge fields do 
not have to be added by hand to the DTQW dynamics, 
as the connection is already part of the basic definition 
of DTQWs and most DTQWs are by definition locally 
gauge invariant m- Second, the difference operators 
(discrete derivatives) which arise in conjunction with the 
local gauge invariance of DTQWs are more complicated 
than the usual finite difference operators used in lattice 
gauge theories. The mathematical properties of discrete 
gauge theories based on DTQWs are thus probably very 
different from the mathematical properties of LGTs. 

Finally, DTQWs are useful in a much wider context 
than high energy or condensed matter physics. DTQWs 
are in particular universal building blocks of quantum al¬ 
gorithms [46] and our results therefore have implications 
for quantum information. For example, the exploration 
of graphs by DTQWs could be influenced by creating 
discrete gauge fields on these graphs. Indeed, not only 
do gauge fields influence the transport of single DTQWs, 
but gauge theories provide a novel manner to implement 
interaction between DTQWs. 
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